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1.1 Introduction 





A comb is a simplified model for various types of natural phenomena which belong to the 
loopless graphs category. The comb consists of a backbone along the horizontal axis and 
fingers or teeth along the perpendicular direction (see Figure fOI for a two sided comb). 

Comb like models have been applied to mimic ramified structures as spiny dendrites 
of neuron cells [MI13, IM13] or percolation clusters with dangling bonds [WH86]. We are 
interested in the first example, where a comb structure with one sided teeth of infinite 
length can be used to describe the movement and binding dynamics of particles inside the 
spines of dendrites. These spines are small protrusions from many types of neurons located 
on the surface of a neuronal dendrite. They receive most of the excitatory inputs and their 
physiological role is still unclear although most spines are thought to be key elements in 
neuronal information processing and plasticity [YuslO]. Spines are composed of a head 1 
p,m) and a thin neck 0.1 p,m) attached to the surface of dendrite (see Fig. II.2|) . 
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FIGURE 1.1 Two sided comb 
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FIGURE 1.2 Electron tomogram of a spiny dendrite. Image taken from Internet 
(http://www.cacr.caltech.edu/projects/Idviz/results/levelsets/ ). 


The heads of spines have an active membrane, and as a consequence, they can sustain 
the propagation of an action potential with a rate that depends on the spatial density of 
spines [EB09]. Decreased spine density can result in cognitive disorders, such as autism, 
mental retardation and fragile X syndrome [NSS02]. Diffusion over branched smooth den¬ 
dritic trees is basically determined by classical diffusion and the mean square displacement 
(MSD) along the dendritic axis grows linearly with time. However, inert particles diffus¬ 
ing along dendrites enter spines and remain there, trapped inside the spine head and then 
escape through a narrow neck to continue their diffusion along the dendritic axis. Re¬ 
cent experiments together with numerical simulations have shown that the transport of 
inert particles along spiny dendrites of Purkinje and Pyramidal cells is anomalous with an 
anomalous exponent that depends on the density of spines [SWDA06, SWDAll, DHll]. 
Based on these results, a fractional Nernst-Planck equation and fractional cable equation 
have been proposed for electrodiffusion of ions in spiny dendrites [HLW08]. Whereas many 
studies have been focused to the coupling between spines and dendrites, they are either 
phenomenological cable theories [HLW08, CBOSb] or microscopic models for a single spine 
and parent dendrite [BS05, SMSOl]. More recently a mesoscopic non-Markovian model for 
spines-dendrite interaction and an extension including reactions in spines and variable resi¬ 
dence time have been developed [FM08, MFHIO]. These models predict anomalous diffusion 
along the dendrite in agreement with the experiments but are not able to relate how the 
anomalous exponent depends on the density of spines [SWDAll, DHll]. Since these exper¬ 
iments have been performed with inert particles (i.e., there are not reaction inside spines 
or dendrites) one concludes that the observed anomalous diffusion is due exclusively to 
the geometric structure of the spiny dendrite. Recent studies on the transport of particles 
inside spiny dendrites indicate the strong relation between the geometrical structure and 
anomalous transport exponents [SWDAll, BBZ09, BWKllj. Therefore, elaboration such 
an analytic model that establishes this relation can be helpful for further understanding 
transport properties in spiny dendrites. The real distribution of spines along the dendrite, 
their size and shapes are completely random [NSS02], and inside spines the spine necks 
act as a transport barrier [BS05]. For these reasons one may reasonably assume that the 
diffusion inside spine is anomalous. In this chapter, we describe some models, based on 
comb-like structure that mimic a spiny dendrite; where the backbone is the dendrite and 
the teeth (lateral branches) are the spines. The models predict: i) anomalous transport 
inside spiny dendrites, in agreement with the experimental results of Ref. [SWDA06], ii) 
also explain the dependence between the mean square displacement and the density of 
spines observed in [SWDAll] and hi) the mechanism of translocation wave of CaMKH 
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(Ca^'*’ - calmodulin-dependent protein kinase II, a key regulator of the synaptic function). 
The chapter is organized as follows. First we study the statistical properties of combs and 
explain how to reduce the effect of teeth on the movement along the backbone as a wait¬ 
ing time distribution between consecutive jumps. Second, we justify an employment of a 
comb-like structure as a paradigm for further exploration of a spiny dendrite. In particular, 
we show how a comb-like structure can sustain the phenomenon of the anomalous diffu¬ 
sion, reaction-diffusion and Levy walks. Finally, we illustrate how the same models can be 
also useful to deal with the mechanism of ta translocation wave / translocation waves of 
CaMKII and its propagation failure. We also present a brief introduction to the fractional 
integro-differentiation in appendix at the end of the chapter. 

1.2 Random walks in combs 


The statistical properties of comb-like structures have been widely studied in the last cen¬ 
tury. The first passage-time and the survival probability were studied by some authors 
[WH86, RedOl]. More recently the interest has been centered in the waiting time distribution 
equivalent to perform a random walk along the teeth [VdB89, CMOS, DBMZ07, MFHIO], 
the mean encounter time between two random walkers [ABC 14] and the occupation time 
statistics [RBI3]. In this section we illustrate two methods to compute the waiting time 
distribution that mimics the effect of a random walk along a teeth. 

We first consider the case of a discrete ID chain where the nearest neighbors are sepa¬ 
rated by a distance a. A random walk, where each walker moves only to one of its nearest 
neighbors with equal probability after a fixed waiting time r, is characterized by the fol¬ 
lowing probability distribution functions (pdf)s for the waiting times and jump lengths, 
respectively, 


(j){t) = 6{t-T), 

w{x) = - [(5(a; — a)-I-<5(a;-I-a)] . (1-1) 

In this way, systems with discrete time and space can be analyzed in terms of the 
Continuous-Time-Random -Walk (CTRW). Next, we add to every site of the backbone a 
secondary branch of length I, to produce a one sided comb-like structure (see Fig. [E3D. 
On such a structure, a walker that is at a given site of the backbone, can spend a certain 
amount of time in the secondary branch before jumping to one of the nearest neighbor sites 
on the backbone. If we are only interested in the behavior of the system in the direction of 
the backbone, then the secondary branches introduce a delay time for jumps between the 
neighboring sites on the backbone. The random walk on the comb structure can be modelled 
as a CTRW with (11.11) and a renormalized waiting pdf (j>(t) that includes the effect of the 
delay due the motion along teeth. 

To determine analytically the effect of the teeth, we invoke convolution rules that were 
introduced in [VdB89] for the case of homogeneous lattices: 

(i) Consider a walker that is initially at a certain site within a tooth. If the walker 
proceeds further into a tooth, i.e., moves away from the backbone, its probability to return 
to the initial site after a time t is a convolution of factors, i.e., a product in the Laplace 
space. 

(ii) The total probability for that walker to return to the initial site is determined by 
summing over all t from 0 to oo. 

(iii) When the walker reaches a crossing, where it can choose between different directions, 
the total probability is the sum of the probabilities for each possible direction. 
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FIGURE 1.3 Sketch of a comb structure with a nearest-neighbor distance a and teeth of length 1 . The 
symbols A, A', A"... denote sites on the lattice (see text). 


The teeth have no other additional crossings. It is as shown in Fig. 11.31 For the sake 
of generality, we consider the case that when the walker is at a site on the backbone, it 
can jump to another site on the backbone with probability ct, or move onto the secondary 
branch with probability 1 — a. 

Without loss of generality, we assume that initially the walker is located on the backbone, 
and we apply the rules (i) - (iii) to determine (/)(t). We study three specific cases. 

(a) Comh structure with I = a. In this case there is only one site on the tooth. A'. The 
walker can only jump in the direction of the backbone with probability a or move onto the 
branch with probability 1 — a and then return to the initial site at the next jump. The time 
it takes to reach one of the nearest neighbors on the backbone is t = t with probability a, 
t = 3t with probability (1 — a) x 1 x a, t = 5t with probability (1 — a)^ x x a, and so 
on. We can write intuitively the general form as 


(l>{t) ay (2j - 1 )t] . (1.2) 

1=1 


The rules listed above for (j){t) should reproduce this behavior. For this purpose we need 
to work in the Laplace space. Let 4>{s) be the Laplace transform of (j){t). The rules (i) - (iii) 
lead to the expression 


OO 

^(s) = a^o [(1 - a)yi 


3=0 


1 - (1 - 


(1.3) 


where (j)o is the probability distribution for a single jump, (j)o = e which is the Laplace 
transform of dnj. 

Equation (11.311 is derived as follows. The term (1 —a)(^g in the sum represents, according 
to rule (i), the probability function for each occurrence of the walker moving onto the 
secondary branch. This expression must be summed up to infinity, according rule (ii), to 
take into account that the walker can move onto the tooth 1, 2,..., oo times. The factor a(f>o 
accounts for the final jump to the nearest neighbor on the backbone. 

It is easy to see that the expression (11.31) may be written as a Taylor series, 

2^* 1 

yis) = {^o) , ( 1 - 4 ) 

1=1 

which is the Laplace transform of dm. The method for determining (j){s) is shown to be 
valid in this case. 
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(b) Comb structure with I = 2a. The secondary branch is two-sites long, A! and AC. 
Similarly to the previous case we can write the distribution for the time probabilities as 


0(s) = a<^o X! 


(1 - a) 


iE 

/c=0 


a(^o(2 - 

2 - (2 - a)^l ■ 


(1.5) 


In this equation a new sum over the index k appears, because the walker can move away 
from the backbone twice. For each such occurrence we must apply rule (i). We also assume 
that jumps to the nearest neighbor occur with probability 1/2 on the linear teeth. 

(c) Comb structure with I ^ oo. Each time the walker moves away from the backbone, a 
new convolution factor appears in 4>{s). For the case I —>■ oo, we have in principle infinitely 


many convolution factors in the expression for 4>{s). Fortunately, we can simplify this situa¬ 
tion considerably. Assume that the walker is at the first site on the secondary branch, point 
AA in Fig. 11.31 and moves away from the backbone. Let r\A' be the probability distribution 
of returning for the first time to the point AA after a time t. Now suppose the same situation 
but for the initial point A!'. It is easy to see that as I ^ oo the limit rjA" —^ tja' has to 
hold, and we can again use the rules (i) - (iii) to determine rjA’- Doing so, we obtain the 
expression 



( 1 . 6 ) 


This expression is equivalent to the form (O with a = 1/2; on the secondary branch 
every jump to a nearest neighbor occurs with probability 1/2. Introducing the condition 
VA" = VA', which is strictly correct for Z = oo, and solving (HH), we find 


VA' = 



(1.7) 


With this result, the distribution 4>{s) is obtained straightforwardly from the three rules (i) 

- (iii): 

rvi . ') 

1J a0o 


^(s) = a^o E [(^ “ a)hVA' 

3=0 


a+{1-a) Jl-(j)l 


( 1 . 8 ) 


In finding the waiting time pdf associated with the comb structures, it is assumed that 
we are only interested in the dynamical behavior along the backbone and the rest of the 
structure is considered as secondary. 

The second method consists in finding the master equation for the random walk moving 
along the backbone. This master equation has to have incorporated the movement along 
the teeth. Consider first the simplest case of a one sided comb with a single node as shown 
in Figure [L4] 

Let a be the probability of moving along the backbone, when the walker is at a node of 
the backbone and 1 — a is the probability of jumping to the teeth, if the walker is at the 
backbone. So, if the movement along the backbone is considered isotropic, the probability 
of jumping to the right or to the left is a/2. The master equation for the probability of 
finding a walker located at the node i of the backbone at time t is (see Figure HTT)) 


t) = ^P{i + a,t-T) + ^P{i - a,t-T)+ P{1, t - t) , (1.9) 

where we have taken into account that the walker waits a constant time r at every node 
between consecutive jumps. The master equation for the node 1 of the tooth reads 


Pi^C) = (1 - a)P{i,t- t) 


( 1 . 10 ) 
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FIGURE 1.4 Comb with a single node along the tooth. In the picture are written the probabilities of 
jumping between nodes of length a. 


that is coupled to (HH). Since time is a discrete variable it is convenient to transform the 
time coordinate to new variable z through the transformation 

OO 

P{z,z) = Y,z*P{z,t). 

Multiplying (HU and (11.101) by z* and summing from t = 0 to infinity we have, respectively 
Pih z) = + a,z) + + bP{l, z) , (1.11) 

P(l,z) = 6(l-ct)P(t,z), (1.12) 

where b = z'^. Inserting (I1.12p into (11.111) and rearranging terms we get the following master 
equation for the movement along the backbone 

ip(f + a,z) + ip(i -a,z) = P{i, z )-—. (1.13) 

It is convenient to stress that (11.131) is actually the master equation for a walker moving 
on a comblike constructed by repeating the element depicted in Figure [TU Finally, we can 
write Eq. (11.131) in the Laplace space for time by taking into account that z = e“®, so that 
b = Hence, (11.131) becomes 

1 « 1 - 1 - 11 - 

-P{z + a,s) + -P(* - a, s) = Pi^, s) --. (1.14) 

2 2 ae 

The example consists in generalizing the above structure to a tooth with R nodes or 
length I = aR. Then, we consider a one sided comb with whose basic element is given in 
Figure [TU The master equation for the probability of finding a walker at the node i of the 
backbone at time t is, at the z space 

z) = Y-P(* + a,z) + ^P{i -a,z) + ^P(l, z). (1.15) 

At the node R the master equation is 

P(P,z) = ^P(P-l,z), (1.16) 

and at the nodes P — 1 and 1 they are 

P{R - 1, z) = ^P(P - 2, z) + 6P(P, z), 


(1.17) 
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FIGURE 1.5 Comb with i? nodes along the teeth. In the picture are written the probabilities of jumping 
between nodes of length a. 


P(l, z) = b(l - a)P(i, z) + ^P(2, z). (1.18) 

At a generic node j of the tooth we get 

P{j, z) = ^P{j -l,z) + ^P{j + l,z), (1.19) 

where j = 2,...,i? — 2. We may solve (11.1911 by proposing the solution P{j,z) = AX^. 
Inserting this solution we obtain the characteristic equation — 2X/b +1 = 0. Hence, 

-f’(+^) = ^1-^++ ^2^'^ , (1-20) 


where A 1.2 are constant to be determined and 

A± = - (1 ± Vl - b^^ . 

On setting j = 1 and j = 2 into Eq. (11.2011 . we get the system of algebraic equations for the 
constants Ai _2 

AiA-i- + A 2 A_ = H(l, z ), 

AiXl+A2Xi=P{2,z). ( 1 . 21 ) 

On the other hand, by inserting (ll.20p into (jl.1611 and (11.1711 we get the equation relating 
Ai and A 2 


Ai 


vfi-2 


3b^ 


1-^ l-TT A 


62 




+A 2 


\R-2 


36= 


1-^ 1-^ A 


62 


vfl-3 


= 0 . 


( 1 . 22 ) 

By combining (jl.22ll and (ll.2ip we can express P{2,z) in terms of P{l,z) in the form 
P(2,z) = h{b)P{l, z), where 


6 ( 6 ) = 


1 

6 


62(A« + A(*) 

A5 + A« + (A« - A5)vT^^ 


2 - 


(1.23) 
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Now, we are in position to get the master equation for the motion along the backbone by 
substituting P(2, z) = h{b)P{l, z) into (11.181) and the result into (11.151) . The final result is 


^P{i + a,z) + ^P{i 




6^(1 — a) 
2-bh{b) 


Pii,z) . 


(1.24) 


We show now how to reduce the effect of the teeth on a waiting time distribution for the 
motion of a walker along the backbone. To this end we make use of the CTRW, and in 
particular, of the generalized master equation for finding a walker at point x at time t when 
it moves a ID space 




J — J P{x — x',t — t')w{x')dx'— P{x,t') 


dt'. 


(1.25) 


Considering the jump lengths distribution given in (ini, for jumps between the consecutive 
nodes of a ID lattice with spacing a and transforming into the Laplace space in time, Eq. 
(jl.25p becomes 


^P(i + a,s) + ^P(i-a,s) 


P(z,s) 


s + K(s) 

^(5) 


(1.26) 


where the memory kernel K{t) is related to the waiting time distribution through their 
Laplace transforms 

K{s) = -^^ (1.27) 

1 - (/)(s) 

Therefore, Eq. (11.261) reduces to 


ip(i + a,s) + ip(f - a, s) = P(f, s)^ 


(1.28) 


The waiting time distribution is obtained by comparing Eqs. (I1.28|) and (I1.24L to get 


ks) 


ab 




1 + 


\R_ \R __■ 


(1.29) 


where b = e If the one sided comb has a teeth with one node, P = 1, thus from Eq. 
(11.291) . the waiting time distribution is 


ks) 


ab 

1 - (1 - a )62 ’ 


which coincides with the result obtained in Eq. (II.3|) . When P = 2, then 


ks) 


ab{2 — 6^) 
2- (2-a)62 ’ 


which is the same result obtained in (Hm. Finally, let us consider the infinite length of the 
teeth P —7> oo, or I a. In this case, Eq. (11.291) reduces to 


ks) 


ab 

Of -I- (1 — a) a/1 — 52 ’ 


(1.30) 


that is equal to (11.81) . In the limit of large times, s —>■ 0, the waiting time distribution in 
Eq. (11.301) reads (j){s) ~ 1 — (1 — a)\/2sT j a, which predicts anomalous diffusion along the 
backbone. So that, the teeth have to have an infinite length to predict anomalous diffusion 
along the whole comb at the asymptotically large times. 















































Comb models for transport along spiny dendrites 


1-9 


smooth dendrite 



FIGURE 1.6 Draw of a smooth dendrite and a spiny dendrites 

1.3 Comb-like models mimic spiny dendrites 


As shown in previous studies, the geometric nature of spiny dendrites plays essential role 
in kinetics [SWDA06, SWDAll, DHll, BKH07, BBZ09, BWKll], The real distribution of 
spines along the dendrite, their size and shapes are completely random [NSS02], and inside 
spines, not only the spine necks but the spine itself acts as a transport barrier [SWDAll, 
BWKll, BS05]. Therefore a reasonable assumption is a consideration of anomalous diffusion 
along both the spines and dendrite. So, we propose models based on a comb-like structure 
that mimics a spiny dendrite, where the backbone is the dendrite and teeth (lateral branches) 
are the spines, see Fig. 11.61 fWe distinguish between a smooth and a spiny dendrite). In this 
case dynamics inside teeth corresponds to spines, while the backbone describes diffusion 
along dendrites. Note that the comb model is an analogue of a ID medium where fractional 
diffusion has been observed and explained in the framework of the CTRW [WH86, AB91, 
MKOO, LZ98] and making use of macroscopic descriptions [ZPRC08] . 

Before embarking for the CTRW consideration in the framework of the comb model, let 
us explain how anomalous diffusion in the comb model relates to the CaMKII transport 
along the spiny dendrite, and how geometry of the latter relates to the anomalous transport. 
As admitted above, the spine cavities behave as traps for the contaminant transport. As 
follows from a general consideration of a Markov process inside a finite region, the pdf of 
lifetimes inside the cavity with a finite volume and arbitrary form decays exponentially 
with time t (see, for example, [RedOl]) (p(t) = ^ exp(—^). Here r is a survival time (mean 
life time), defined by the minimum eigenvalue of the Laplace operator and determined by 
geometry of the cavity. For example, in Refs. [BBZ09, BWKll], for spines with a head of 
volume V and the cylindrical spine neck of the length L and radius o, the mean life time is 
r = LV/ira^D = LF'/D, where D is diffusivity of the spine. Therefore, the mean probability 
to find a particle inside the spine after time t, (i.e, the survival probability inside the cavity 
from 0 to t) averaged over all possible realizations of r is given by the integral 

pOO pOO 

= / / ’fit'/T)f{T)dTdt' (1.31) 

Jt Jo 

where /(t) is a distribution function of the survival times r (recall that size and shape of 
spines are random [NSS02]). Finally, the waiting time pdf can be easily calculated from Eq. 
(I1.3ip . as follows 

poo 

(j){t) = -dt'i’it) = / (p{t/T)f{T)dT. (1.32) 

Jo 

In the simplest case, when the distribution is the exponential /(r) = (I/tq) exp(— t/tq). 
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one obtains from Eq. (I1.32|) that the general kinetics is not Markovian and the waiting time 
pdf is a stretched exponential for large times 


(/>(t) = 


-f 

To Jo 


OO t/T 


= ^Kq (2^/t/To^ ^ exp(-v'V'^o) ,t/To » 1 . 


The situation is more interesting, when the distribution of the survival times is the power 
law /(r) ^ (0 < a < 1). In this case the waiting time pdf is the power law as 

well (j){t) ^ that leads to subdiffusion motion along the dendrite. This result follows 

from the CTRW theory, since all underlying micro-processes are independent Markovian 
ones with the same distributions [MKOO]. 

Now we explain the physical reason of the possible power law distribution 4>{t). At this 
point we paraphrase some arguments from Ref. [bAH05] with the corresponding adaptation 
to the present analysis. Let us consider the escape from a spine cavity from a potential 
point of view, where geometrical parameters of the cavity can be related to a potential U. 
For example, for the simplest case, mentioned above, it is t/ = VL/tto^, which “keeps” 
a particle inside the cavity, while Dtq plays a role of the kinetic energy, or the “Boltz¬ 
mann temperature”. Therefore, escape probability from the spine cavity/well is described 
by the Boltzmann factor e:xjp{—U/ Dtq) . This value is proportional to the inverse waiting, 
or survival time 

As admitted above, potential U is random and distributed by the Poisson distribution 
P{U) = Uq^ e^p{—U/Uo), where Uo is an averaged geometrical spine characteristic. The 
probability to find the waiting time in the interval (t, t + dt) is equal to the probability to 
find the trapping potential in the interval {U, U + dU), namely cj){t)dt = P{U)dU. Therefore, 
from Eq. (11.331) one obtains 

~ ^ • (1-34) 

Here 7 = € (0,1) establishes a relation between geometry of the dendrite spines and 

subdiffusion observed in [SWDA06, SWDAll] and support application of the comb model, 
which is a convenient implement for analytical exploration of anomalous transport in spiny 
dendrites in the framework of the CTRW consideration. 


1.3.1 Anomalous diffusion in spines 

Geometry of the comb structure makes it possible to describe anomalous diffusion in spiny 
dendrites structure in the framework of the comb model. 

Usually, anomalous diffusion on the comb is described by the 2D distribution function 
P = P{x, y, t), and a special behavior is that the displacement in the x-direction is possible 
only along the structure backbone (a;-axis at y = 0). Therefore, diffusion in the x-direction is 
highly inhomogeneous. Namely, the diffusion coefficient is D^x = Dx6(jj), while the diffusion 
coefficient in the y-direction (along teeth) is a constant Dyy = Dy. Due to this geometrical 
construction, the flux of particles along the dendrite is 

Jx = -DxS{y)^ (1.35) 


and the flux along the finger describes the anomalous trapping process that occurs inside 
the spine 




dP 

RL 


(1.36) 
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where P{x,y,t) is the density of particles and 


fit) = p^fit) (1.37) 

RL 

is the Riemann-Liouville fractional derivative, where the fractional integration If is defined 
by means of the Laplace transform 

C[Iffit)] = s-'>fis). (1.38) 


01-7 

dt^'y 


So, inside the spine, the transport process is anomalous and (y^(t)) C, where 7 G (0,1). 
Making use of the continuity equation for the total number of particles 

^-tdivJ = 0, (1.39) 


where J = (J^,, Jy), one has the following evolution equation for transport along the spiny 
dendrite 


dP 

dt 


DJiy) 


d^P 

dx'^ 


-A, 


91-7 

dP-^ 


d^P 

RL dy'^ 


= 0 . 


(1.40) 


The Riemann-Liouville fractional derivative in Eq. (I1.40|) is not convenient for the Laplace 
transform. To ensure feasibility of the Laplace transform, which is a strong machinery for 
treating fractional equations, one reformulates the problem in a form suitable for the Laplace 
transform application. 

To shed light on this situation, let us consider a comb in the 311 [bAH05]. This model 
is described by the distribution function Piix,y, z,t) with evolution equation given by the 
equation 


dPi 


DJ{y)5{z) 


d^Pi 

dx"^ 


Dy5{z) 


d'^Pi 

W 


d^Pi 


= 0. 


(1.41) 


It should be stressed that z coordinate (do not confuse with the z variable introduced 
in the previous section) is a supplementary, virtue variable, introduced to described frac¬ 
tional motion in spines by means of the Markovian process. Thus the true distribution is 
P{x, y, t) = Pi (x, y, z, t)dz with corresponding evolution equation 


dP 


DxSiy) 


d^Piiz = 0) ^ d^Piiz = 0) 


dx^ 


-A 


9y2 


= 0 . 


(1.42) 


A relation between P(x, y, t) and Pi(x, y, z = 0, t) can be expressed through their Laplace 
transforms 

Pi(x,y,z = 0,s) = ^P(x,y,s), (1.43) 

where P{x,y, s) = £[P(x, y,t)] and Pi{x,y,z,s) = £[Pi{x,y, z,t)]. 

Equation (11.4311 establishes a relationship between the distributions Pi(x,y, z = 0,t) 
and P(x, y, t) in the Laplace space. Both distributions are related through the expression 


/ OO 

Pi ix,y,z,t) dz. 

-00 


If we transform the above equation by Fourier-Laplace we get 

P { kxf ky^ s) -Pi { kxt ky^ kz 0; '5) • 


(1.44) 
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Then, Eq. (11.431) is nothing but a relation between Pi {kx, ky, kz = 0, s) and Pi {kx, ky, z = 0, s). 
To find Pi{kx,ky,kz,s) we transform Eq. (11.411) by Eourier-Laplace and after collecting 
terms we find 


{J^Xtky^kz^ — 


1 - DxklPi {kx,y = 0 ,z = 0 , s) - Dyk^Pi {kx,ky,z = 0 ,s) 

s + k'i 


(1.45) 


where the initial condition has been assumed Pi (x, y, z,t = 0) = d (x) S (y) S (z) for simplic¬ 
ity. Setting kz = 0 one gets 


Pi ( kx , ky ^ kz — 0 ; ^) 


1 - DxklPi {kx, y = 0 , z = 0 , s) 

s 


Inverting Eq. (11.451) by Fourier over kz we obtain 


PykyPi {kx, ky, z — 0 , s) 

(1.46) 


Pi {^X , ky , Z , s') - 


1 - DxklPi {kx,y = 0 , z = 0 , s) - Dyk^Pi {kx,ky,z = 0 , s) 




Then setting z = 0, one obtains 


Pi ( kx, ky, z — O^a) — 


1 - DxklPi {kx, y = 0 , z = 0 , s) - Dyk^Pi {kx,ky, z = 0 , s) 

^s 


(1.47) 


Combining Eqs. (11.461) and (I1.47L one has 


Pi {kx, ky, Z 0, s) 2 Pi {kx, ky, kz 0; 5) , 

then the Fourier inversion over kx and ky yields Eq. (11.431) . Finally, performing the Laplace 
transform of Eq. (11.421) one obtains 


sP{x,y,s) - DxS{y) 


d^Pi{x,y,z = 0,s) d^Pi{x,y,z = 0,s) 


dx"^ 


- D 


5y2 


= P{x,y,t = 0), (1.48) 


and substituting relation (11.431) , dividing by y/s and then performing the Laplace inversion, 
one obtains the comb model with the fractional time derivative 


dip 


32 p g2p 


(1.49) 


where 2Dx^y —>■ Dx^y and the Caputo derivativ^B can be defined by the Laplace transform 
for 7 e (0,1) [Mai96] 


= sy{s) - s'^-^f{t = 0). (1.50) 

The fractional transport takes place in both the dendrite x direction and the spines y co¬ 
ordinate. To make fractional diffusion in dendrite normal, we add the fractional integration 


dP 


*To avoid any confusion between the Riemann-Liouville and the Caputo fractional derivatives, the former 
one stands in the text with an index RL: while the latter fractional derivative is not indexed 

^ 3 -. Note, that it is also convenient to use Eq. 111.501 1 as a definition of the Caputo fractional derivative. 
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ll ^ by means of the Laplace transform (11.381) . as well C 
Eq. (jl.49|) . after generalization ^ 7 G (0,1), 


ri-7 


dfy 


dx"^ dy'^ 


fit) 


= 0 . 




Performing the Fourier-Laplace transform in (11.511) we get 


P ( hx -/kyj s') — 


Pjkx, ky, t = 0)- DxklPjkx, y = 0,s) 
s + Dyk^s^-^ 


/(s). This yields 

(1.51) 

(1.52) 


where the Fourier-Laplace image of the distribution function is defined by its arguments 
CPxiFy[Pix,y^t)] = P{kx,ky,s). If P{x,y,t = 0) = S(x)d{y), inversion by Fourier over y 
gives 


Pikx,y,s) = 


1 - DxklP{kx,y = 0 ,s) 


Taking y = 0 the above equation provides 

P{kx,y = 0 , s) = 


exp 




s(2-7)/2y;D;+Zl,fc2 ’ 


which yields after inserting into Eq. (11.521) 

Pikx: ky, s) = ■ 

s + Dyk^s^ ^ 


1 - 


n 

rJxi^x 


S<^^-^)/^^y + Dxkl 


(1.53) 


(1.54) 


(1.55) 


We can calculate the density of particles at a given point x of the dendrite at time t, namely 
P{x,t), by integrating over y in the Fourier space 


then 


so that 


ky - O7 5 ) - 


(1.56) 

5(2-7)72/0;+ ’ 


2Dx 1 

(1.57) 

/c,=o 


(1.58) 


Equation (11.581) predicts subdiffusion along the spiny dendrite that is in agreement with 
the experimental results reported in [SWDA06]. It should be noted that this result is coun¬ 
terintuitive. Indeed, subdiffusion in spines, or fingers should lead to the slower subdiffusion 
in dendrites, or backbone with the transport exponent less than in usual comb, since these 
two processes are strongly correlated. But this correlation is broken due to the fractional 
integration ll~'^ in Eq. (11.511) . On the other hand, if we invert (I1.56P by Fourier-Laplace we 
obtain the fractional diffusion equation for P{x,t) 

d^-ip _ Dx d^P 

dC~i dx"^ 

which is equivalent to the generalized Master equation (11.251) . in the diffusion limit: 


dP 




(1.59) 
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if the Laplace transform of the memory kernel is given by K{s) = which corre- 

spends to the waiting time pdf in the Laplace space given by 


<^(s) 



(1.60) 


that is (pit) ^ as t ^ oo. Let us employ the notation for a dynamical exponent dw 

used in [SWDA06, SWDAll]. If dw = 4/(2 — 7 ) then the MSD grows as Qn the other 
hand, it has been found in experiments that du, increases with the density of spines ps and 
the simulations prove that d^ grows linearly with pa. Indeed, the experimental data admits 
almost any growing dependence of d^ with p^ due to the high variance of the data (see Fig 
5.D in [SWDAll]). Equation (|1.58l) also establishes a phenomenological relation between 
the second moment and pa. When the density spines is zero then 7 = 0, du, = 2 and normal 
diffusion takes place. If the spine density ps increases, the anomalous exponent of the pdf 
(11.601) 1 — 7/2 = 2/dw must decrease (i.e., the transport is more subdiffusive due to the 
increase of ps) so that dyj has to increase as well. So, our model predicts qualitatively that 
dw increases with ps, in agreement with the experimental results in [SWDAll]. 


1.3.2 Levy walks on fractal comb 

In this section we consider a fractal comb model [lomll] to take into account the inhomo¬ 
geneity of the spines distribution. Here, we consider the comb model for a phenomenological 
explanation of an experimental situation, where we introduce a control parameter which 
establishes a relation between diffusion along dendrites and the density of spines. Suggesting 
a more sophisticated relation between the dynamical exponent and the spine density, we 
can reasonably suppose that the fractal dimension, due to the box counting of the spine 
necks, is not integer: it is embedded in the ID space, thus the spine fractal dimension is 
V S (0,1). According the fractal geometry (roughly speaking), the most convenient param¬ 
eter is the fractal dimension of the spine volume (mass) /J.spine(a^) = pix) ^ \x\''■ Therefore, 
following Nigmatulin’s idea on a construction of a “memory kernel” on a Cantor set in the 
Fourier space [Nig92] (and further developing in Refs. [LMNN, RLWQ03, Bill]), this 

leads to a convolution integral between the non-local density of spines and the probabil¬ 
ity distribution function Pix,y,t) that can be expressed by means of the inverse Fourier 
transform [lomll] [\kx\^~''P{kxTy,t)]. Therefore, the starting mathematical point of 
the phenomenological consideration is the fractal comb model 

Dx5iy)I^ \\kx\^ ''P{kx,y.,t)\ =Q . (1-61) 

Performing the same analysis in the Fourier-Laplace space, presented in previous section, 
then Eq. (11.5611 reads 


P{kx, s) = P{kx,ky = 0 , s) 


si-^-PP^y + Dx\kx\P ’ 


(1.62) 


where jd = 3/2 -|- v 

Contrary to the previous analysis expression (11.5711 does not work any more, since su- 
perlinear motion is involved in the fractional kinetics. This leads to divergence of the second 
moment due to the Levy flights. The latter are described by the distribution ~ , 

which is separated from the waiting time probability distribution (pit). To overcome this de¬ 
ficiency, we follow the analysis of the Levy walks suggested in [ZKB90, ZK95]. We consider 
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our exact result in Eq. (11.621) as an approximation obtained from the joint distribution of 
the waiting times and the Levy walks. Therefore, a cutoff of the Levy flights is expected 
at |a;| = t. This means that a particle moves at a constant velocity inside dendrites not all 
times, and this laminar motion is interrupted by localization inside spines distributed in 
space by the power law. 

Performing the inverse Laplace transform, we obtain solution in the form of the Mittag- 
Leffler function [BH55] 

P(k,,t) = , (1.63) 

where D = —Ss=. For the asymptotic behavior |fc| —>■ 0 the argument of the Mittag-Leffler 

V 

function can be small. Note that in the vicinity of the cutoff |x| = t this corresponds to the 
large t (|fc| ^ \ <C 1), thus we have [BH55] 


Ei -^/2 [-D\k\^t^ Ri exp 


D\k\PC-^/^\ 
r(2-7/2); ■ 


Therefore, the inverse Fourier transform yields 


(1.64) 


Pix,t) ^ A,-p(2_^/2)|^|(5+.)/2 ’ 


(1.65) 


where is determined from the normalization condition (The physical plausibility of 
estimations (jl.64|) and (jl.65l) also follows from the plausible finite result of Eq. (11.6511 . 
which is the normalized distribution P{x,t) where \x\ = t). Now the 

second moment corresponds to integration with the cutoff at x = t that yields 

{x^{t)) = , ( 1 . 66 ) 

where i, = -^ ^ generalized diffusion coefficient. Transition to the 

absence of spines means first transition to normal diffusion in teeth with 7 = 1 and then 
1 / = 0 that yields 

{x\t)) = K^^ot. (1.67) 


1.3.3 Fractional reaction-diffusion along spiny dendrites 


Geometrically, spiny dendrites in the HD space are completely described by a comb structure 
in the 2D, where the spine density on the cylinder surface is projected on the II? axis (say 
the X axis): p(x,r = const, 6 () —)> p(x). Here p(x,r = const,9) is the spine density, while 
p{x) is the density of the comb teeth. In what follows, we consider p(x) = g = const, which 
is, probably, the most realistic case. Fractional diffusion inside the spines is described by 
fractional diffusion inside the teeth. Therefore, one considers a two-sided comb model as in 
Figure 11.11 and the starting mathematical point of the phenomenological consideration is 
the Fokker-Planck equation obtained in [Mil 3]. 

It reads. 


dP 


32 p Q2p 


dx'^ 


dy^ 


( 1 . 68 ) 


This equation is obtained by the re-scaling with relevant combinations of the comb pa¬ 
rameters Dx and Dy, such that the dimensionless time and coordinates are D^t/Dy — t 
Dx/Dy —>■ X, Dx/Dy — >■ y/y/g, correspondingly [IB05], and parameter g can be considered 
as a constant density of the fingers. 
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As admitted above a variety of interactions inside spines leads to correlated noises in 
dendritic spines [ZHIO]. The strong correlations of that leads to anomalous (subdiffusive) 
motion inside the spines. Following a phenomenological description by the CTRW, this 
subdiffusion is controlled by a waiting-time pdf (^(t) decaying according to the power law. 
Therefore, normal diffusion of the contaminant density P{x,y^t), for example activated 
CaMKII, in spines is replaced by the anomalous transition term 

where K{t) is again the time memory kernel of the generalized master equation and is 
defined in (11.271) . For subdiffusion, 0(<) = with 0 < 7 < 1 that yields [FSMSS08] 

K{s) = 

One may recognize that Eq. (I1.69|) is a formal expression for the anomalous transport 
with very complicated form in the time domain, which in turn, is very inconvenient for the 
analytical treatment. Therefore, the comb model may be presented in the Laplace domain. 
Substituting Eq. (I1.69|) in Eq. (11.6811 . then performing the Laplace transform and taking 
into account Eq. (11.271) . one obtains the comb model in the Laplace domain 


sP = S{y) 


92 P 


gs 


l-a 


d^P 

W 


Pn 


(1.70) 


Here Pq = P{x,y,t = 0) is the initial condition. As admitted, the kernel K{s) is problem¬ 
atic for the Laplace inversion, since it leads to the appearance of the initial condition. To 
overcome this obstacle, one multiplies Eq. (11.701) by and then perform the Laplace 
inversion that yields 


^0 


dP{x,y,t') 

dt' 


- 5{y) 


d^P{x,y,t') 

dx'^ 


dt' = g 


d'^P{x,y,t) 


(1.71) 


Amending this equation by reaction term, one arrives at the integro-differential equation: 




dP{x,y,t') 

dt' 


- Siy) 


d'^P{x,y,t') 


dx"^ 


dt' = g 


d'^P{x,y,t) 

dy'^ 


+ gC[Pix,y,t)] , (1.72) 


which describes the 2D inhomogeneous reaction diffusion in the dispersive medium. Here 
C[P{x,y,t)] = C{P) is a reaction kinetic term. In particular, to model reaction kinet¬ 
ics inside dendrites, it can be considered either linear C{P) = CP, or logistic C'(P) = 
CP{x, y, t)[l — P{x, y, t)] [BrelS]. Integration with the power law kernel D'' ensures anoma¬ 
lous diffusion in both the dendrite and spines. 

In what follows we use convenient notations of fractional integro-differentiation given 
in Eq. (11.501) and the text below. Owing to this notation, from Eq. p.72D the equation for 
P{x,y,t) reads 

= (1.73) 

To give a first and brief insight on the problem of the front propagation, let us consider 
the linear reaction and 7 = 1. In this case, one obtains a “simple” solution [Ioml2] for 
the travelling wave along the x axis (inside dendrites). Introducing the total probability 
distribution function Pi{x,t) = f dyP(x,y,t), one obtains 


Pi{x,t) = exp 

^ ^ TTV^ 

This yields the coordinates of the front x ^ t^P that spreads with the decaying velocity 
V ~ t~^P. This solution illustrates the asymptotic failure of the reaction-transport front 
propagation due to subdiffusion inside spiny dendrites. 


_2y/^ 


-Cgt 


(1.74) 
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1.4 Front propagation in combs 


Recently, a mechanism of translocation wave of CaMKII has been suggested [EBIO]. As 
shown, activated CaMKII contaminant travels along dendrites with additional translocation 
inside spines. Process of activation (the conversion of primed CaMKII to its active state) 
corresponds to the irreversible reaction that, in absence of spines, is described by the Fisher- 
Kolmogorov-Petrovskii-Piskunov (FKPP) equation (It relates also to the logistic kinetic 
term) [FBIO]. Therefore, in the framework of the suggested above scheme of the dispersive 
subdiffusive comb (I1.73L nonlinear reaction at subdiffusion in dendrites takes place along 
the a;-axis bound, while subdiffusion in fingers describes the translocation inside spines. 
Therefore, reaction-transport equation (11.731) now reads 


gip 

dfy 


mii 


1-7 


d‘^P 


- 


(1.75) 


Here D describes the diffusivity inside dendrites, while C{P) = CP{1 — P) is the nonlin¬ 
ear reaction term. Again, integrating over y to obtain the kinetic equation for the total 
distribution Pi(a:,t), we have 


dC ‘ 


1-7 




= 0 . 


(1.76) 


For the brevity, we denoted Pq = P{x, y = 0, t). Consider the fractional comb model (II.75|) 
without reaction, 

dc ^ dx^ ^dy^^' 

and perform the Laplace transform, thus one obtains 

Q2p 


(1.77) 


s'<P-5{y)Ds^- 


dx'^ 


= + ^<5(x)(i(y). 


(1.78) 


where for the initial condition we take P{t = 0) = S{x)6{y). Looking for the solution in 
form 

P{x,y,s) = exp[-\/sT'/y|y|]/(a;,s). (1.79) 

one can see that Pq = P{x,y = 0,s) = f{x,s) and integration Fq. p.79l) over y yields 


Po(x,s) = W —Pi(x,s). 

V 4y 

Substituting (11.801) into (11.761) in the Laplace space, one gets 


(1.80) 


= 0) -110 - 0A = -cc[pSi. (1,81) 

Multiplying this equation by e®* and using identity e®‘s“/(s) = ^e®*s““^/(s), we integrate 
with the corresponding contour to obtain the inverse Laplace transform. This yields 


dPi 


dt 2y7 dt 


1 9 ji-i 




0 

4y 


d 1-2 


The nonlinear term is obtained by the following chain of transformations 


C[P 2 ] = C 


P-'Pn 


= ^ \c-hiPi 

4g I 


7 _ C 


d 1-2 


1 2 


(1.82) 


(1.83) 
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Note that a specific property of these transforniations is an irreversibility with respect to 
the Laplace transform, since, as well known, the Laplace transform of the Riemann-Liouville 
fractional derivative involves the (quasi) initial value terms like Piit = 0) = 5{x) [MKOO]. 

To evaluate the overall velocity of the asymptotic front, let us introduce a small param¬ 
eter, say e, at the derivatives with respect to the time and space [Fre96]. To this end we 
re-scale x —>■ x/e and t —>■ t/e, and Pi{x,t) —7> Pf{x,t) = Pi (|, ^). Therefore, one looks for 
the asymptotic solution in a form of the Green approximation, 


Pi{x, t) = exp 


G^{x,t) 

e 


(1.84) 


The main strategy of implication of this construction is the limit e —>■ 0, one has 
= 0, except the condition when G^{x,t) = 0. This equation determines 


G^{x,t) 


exp 

the position of the reaction spreading front, see Eq. (11.741) . Moreover, we consider the limit 
G{x,t) = lim£_>.o G'^(x, t) as the principal Hamiltonian function [Fre96] that makes it pos¬ 
sible to apply the Hamiltonian approach for calculation of the propagation front velocity. 
In this case partial derivatives of G(a:, t) with respect to time and coordinate have physical 
senses of the Hamiltonian and momentum: 


dG{x, t) 
Wt 


= —P[ and 


dG{x, t) 
dx 


= p. 


(1.85) 


Now the method of the hyperbolic scaling, explained above, can be applied. Therefore, 
we have the ansatz (jl.84l) for the probability distribution function inside dendrites. Inserting 
expression (jl.84l) in Eq. (11.8211 . one considers fractional integrations in time. Let us start 
from the last term in Eq. (I1.82L which is the reaction term. We rewrite it in the following 
convenient form 


d 


r(i-i)9t Jo 

Then performing expansion 


exp[—G®(t — et',x — ex')/e]. 


( 1 - 


, i\ \ dG^{x,t) . dG'^(x,t) . 

G {t — et, X — ex ) ~ G (x, t) — e - —-t — e -- -x 


dt 


dx 


and substituting this in Eq. (I1.86L one obtains 


r(i-i) L 


[ dG^{x,t)) 


r G7a;,t)] 

dt 

exp 

e 


{t') 


/\- 7/2 


exp 


J JO 


dG^{x,t) , dG^x,t) , 

■t H-- X 


at 


dx 


dt'. 


(1.87) 

It should be noted that we neglect differentiation of the upper limit of the integral, since 
this term is of the order of 0(e^'''^/^) ~ o{e) that vanishes in the limit e —>■ 0. The same 
procedure of expansion is performed for the diffusion term in Eq. (11.8211 that yields 


^3 


r(l - ^) dtdx^ 


exp 


G^{x,t) 


(tO’^exp 


dG^{x,t), aG^{x,t) ; 

■t H-- X 


at 


dx 


dt'. (l.S 


Differentiating in the limit e —>■ 0 and taking into account that the Hamiltonian H and the 
momentum p in Eq. (11.851) are independent of x and t explicitly (that leads to the absence 
of mixed derivatives), one obtains the Laplace transform of the subdiffusive kernel t~^. 
After these procedures in Eqs. (11.871) and (11.881) . the kinetic equation (11.821) becomes a kind 
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of Hamilton-Jacobi equation that establishes a relation between the Hamiltonian and the 
momentum 


= 


'Dp^ + C 
_ 2Vff . 


2 

2-y 


(1.89) 


and the action is G{x,t) = /Q[p(s)i:(s) — H{p{s),x{s))]ds. The rate v at which the front 
moves is determined at the condition G{x,t) = 0. Together with the Hamilton equations, 
this yields 

dH H , , 

V = X = , V = — . (1.90) 

op p 

Note that the first equation in (11.901) reflects the dispersion condition, while the second one 
is a result of the asymptotically free particle dynamics, when the action is G(a;, t) = px — Ht. 
Taking into account x = vt, one obtains Eq. (11.901) (see also details of this discussion e.g. 
in Refs. [CFM08, MCF04]). Combination of these two equations can be replaced by 


. H . H{p) 

V = mm ——7 = mm- 

H>0 p(H) p>0 p 


(1.91) 


We also have from the front velocity conditions (11.901) -^\nH = 1/p that, eventually, yields 
from Eq. (|1.9ip 


2 2 + t 

/4\2-^ U / C 
V5/ 2-7 V2-h7j 


(1.92) 


To proceed, we first, admit that the limiting case of this result with 7 = 0 corresponds to 
the CaMKH propagation along the dendrite only (i.e., there are no teeth). Therefore, Eq. 
(11.921) after rescaling D/^ —>• D and G/y/g —>■ G recovers the FKPP scheme for 7 = 0 
that yields v = V DC. 

It should be admitted the absence of the failure of the activation front propagation. It 
has a simple explanation due to the absence of a reaction “sink” term —hP in Eq. (II.75|) 
by neglecting the possibility of spines to accumulate a large amount of Ca^^ [KS06, SegOS], 
where h is a translocation/accumulation rate [EBIO]. Introducing this term in Eq. (11.751) . 
our anticipation is that the hyperbolic scaling for this new equation yields a solution similar 
to Eq. (11.911) with H = 0 that corresponds to the failure of the front propagation. Moreover, 
this asymptotic solution for Pi{x,t) always takes place, as one of possible solutions. 

Inserting the sink term in Eq. (11.751) . one obtains 


dip 

dti 




d'^P 


d'^P 

9-^ - ghP = 0. 


(1.93) 


Repeating the same procedures of the Laplace transform and integration over y with defi¬ 
nition Pi = P(x, y, s)dy, and using the substitute 


Pi {x, s) = 2^g/s"/P{x, y = 0,s), 


one obtains 


sPi — 6{x) 


dPi Cy/7^ - 


hgs^-'^Pi. 


(1.94) 


Here we neglect the nonlinear term, since, as follows from the above analysis, in the further 
hyperbolic scaling approximation, this term does not contribute to the Hamilton-Jacobi 
equation, and we also know how to handle it. Again multiplying this equation by e®* and 
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using the same identity e®‘s“/(s) = as above, we obtain the inverse Laplace 

transform. Thus, Eq. (11.941) reads 


dPi ^ rm c 

dt 2y/g dt ‘ _dx^ D ^ 


hg^nPi 

^ dt ‘ 


Application of the hyperbolic scaling with the asymptotic solution (11.8411 yields 


(1.95) 


2^H 


+ CPi - 2hg5ifi-7 


(1.96) 


Let us consider a specific case 7 = 2/3 that yields 


P = 


3 

+ C - 2hgi'\ " 


(1.97) 


For C > 2hgi there is no failure and the front asymptotically propagates with a constant 
velocity. For C < 2hg'^ the only solution is P = 0 and yields ?; = 0. So, 2hg'^ is the minimum 
reaction rate necessary to sustain propagation along the spiny dendrite due to the presence 
of translocation. Analogously, (( 7 / 2 / 1 )^/^ can also be viewed as the minimum value for the 
density of spines necessary to have propagation failure. Both results are in agreement with 
the results obtained from very different models based on the cable model [CB03a]. 

In general case, one compares the interplay between the activation CP^ and the translo¬ 
cation —2hg^P^~'^ terms in Eq. (11.961) in the limit —>■ 0. For 7 G 1), the translocation 
term is dominant and leads to the solution with P — 0 and the failure of the front propaga¬ 
tion, correspondingly. When 0 < 7 < |, the activation in dendrites can be dominant. This 
situation is more complicated, and the activation-translocation front can propagate with an 
asymptotically finite velocity. 

Finally, let us consider the linear counterpart of Eq. (11.7511 with the linear reaction term 
C(P) = CP. This analysis will be useful to uncover the behavior of the tail of the total 
distribution and check if the front accelerates or travel with constant velocity. Rewrite the 
equation for the total distribution Pi{x,t). As follows from the fractional differentiation of 
Eq. (I1.82|l . this equation reads (see also [MI13]) 


d^-iPi 

dC-i 


D r d^ 
2y/g _dx^ 



(1.98) 


with the initial condition Pi{x,t = 0) = S{x). After the Fourier transform P[Pi{x,t)] = 
Pi{k,t), one obtains the solution in the form of the Mittag-Leffler function 


Pi{k,t) = Ei_x 


A{k)f-i 


(1.99) 


where A{k) = (C — Dk"^)/2^. At the asymptotic condition, when 1, we have 

C ^ Dk^ that yields asymptotic behavior of the Mittag-Leffler function as a growing 
exponent (for the large positive argument) [BH55] 


Pi{k,t) 




After the Fourier inversion, one obtains 


exp 


C A 2 Dk‘^\ 

I 2-7 C ) 


t . 

( 1 . 100 ) 


Pi(x, t) = exp 



2 

2-7 


t 


(2-7)x^(2Ag)^ 

SDC^t 


( 1 . 101 ) 
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that, finally, yields the nonzero and constant overall velocity of the reaction front propaga¬ 
tion. Note that for normal diffusion, 7 = 0, one arrives at the Fisher velocity v = yjDC/g —>■ 
\fDC, see limiting case 7 = 0 in Eq. (11.921) . 


1.5 Conclusion 


In this chapter we show that a comb is a convenient model for analytical exploration of 
anomalous transport and front propagation phenomena along spiny dendrites. We have 
studied the properties of a random walk motion along the backbone in presence of teeth. 
Teeth are lateral branches crossing the backbone and we have shown here that their effect 
on the movement along the whole structure can be reduced to a waiting time distribution at 
the nodes of the backbone during the movement of the random walker. Recent experiments 
and numerical simulations have predicted anomalous diffusion along spiny dendrites. 

We have shown here that this anomalous phenomenon can be explained in the framework 
of a comb model with infinitely long teeth. Moreover, due to the random distribution of 
spines along the parent dendrite and the presence of binding reaction inside spines one can 
present a physically reasonable justification of the power law of the waiting time distribution 
that leads to subdiffusion in both the teeth and the backbone. 

We have shown how to predict anomalous diffusion in spines by constructing a fractional- 
diffusion equation. By using the CTRW formalism we have computed the mean square dis¬ 
placement for the transport along the whole comb. We presented an illustration of how 
to take into account the inhomogeneous distribution of spines along the dendrite by using 
a fractal comb. We have also constructed the corresponding fractional-diffusion equations 
and computed the mean square displacement. From the other hand, the constructed toy 
models are simple enough, like the comb model that makes it possible to suggest and under¬ 
stand a variety of reaction-transport schemes, including anomalous transport, by applying 
a strong machinery of fractional calculus and hyperbolic scaling for asymptotic methods. 
This approach allows to suggest an analytical description of reaction-transport scenarios 
in spiny dendrites, where we consider both a linear reaction in spines, see Eqs. (11.721) and 
(11.731) . and nonlinear reaction along dendrites, considered in the framework of the FKPP 
scheme [Fis37]. To this end we suggest a fractional subdiffusive comb model, where we 
apply a Hamilton-Jacobi approach to estimate the overall velocity of the reaction front 
propagation. We proposed an alternative approach of a recently suggested mechanism of 
translocation wave of CaMKII [EBIO], where activated CaMKII contaminant travels along 
dendrites with additional translocation inside spines, and process of activation corresponds 
to the irreversible reaction described by the FKPP equation (11.821) . One of the main effect, 
observed in the framework of the considered model, is the failure of the front propagation 
due to either the reaction inside spines, or interaction of reaction with spines. In the first 
case the spines are the source of reactions, while in the latter case the spines are a source 
of damping, for example they act as a sink of an activated contaminant (CaMKII). The 
situation is controlled by three parameters CaMKII activation C, CaMKII translocation 
rate h and the fractional transport exponent 7 . The latter reflects the geometrical structure 
of the transport system: when 0 < 7 < |, the activation in dendrites can be dominant, 
and the activation-translocation front can propagate with an asymptotically nonzero and 
constant velocity. For 7 = 2/3 we have found a criteria for the emergence of propagation 
failure or for the sustain of the propagation in terms of the reaction rate, the translocation 
rate and the spine’s density. 

It should be admitted, in conclusion, that physical arguments suggested above, explain 
why anomalous transport, namely subdiffusion, of either CaMKII or neutral particles is 
possible and support implementation of the comb model. These arguments are based on 
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geometry of dendritic spines that determines an expression for the transport exponent in 
Eq. (I1.34p . This situation becomes more sophisticated in a case of the nonlinear FKPP re¬ 
action. Indeed, as shown, the power law kernel of the transition probability considered due 
to the geometry arguments is insensitive to the nonlinear reaction. This consideration dif¬ 
fers completely from a mesoscopic non-Markovian approach, developed in [FM08, MFHIO], 
where spines-dendrite interaction and an extension including reactions in spines have been 
described in framework of variable residence time. This leads to essential complication of 
the transition probability due to the nonlinear reaction term [FedlO, MFHIO]. 

1.6 Appendix: Fractional integro—differentiation 


The consideration of a non-Markovian process in the framework of kinetic equations leads 
to the study of the so-called fractional Fokker-Planck equation, where both time and space 
processes are not local [MKOOj. In this case, derivations are substituted by integrations with 
the power law kernels. One arrives at so-called fractional integro-differentiation. 

A basic introduction to fractional calculus can be found, e.g., in Ref. [Pod98]. Fractional 
integration of the order of a is defined by the operator 


= f(T)it-T)°‘ ^dr, (a>0) 

r(a) Ja 


( 1 . 102 ) 


where r(a) is a gamma function. There is no constraint on the limit a. In our consideration, 
0 = 0 since this is a natural limit for the time. A fractional derivative is defined as an inverse 
operator to alt — ^ 

form is convolution 


= It °' = correspondingly J“ = = Dt°‘. Its explicit 


n“ =_i_ 

‘ r(-a) 


/M 


-dr. 


(1.103) 


/o (i-T)“+i 

For arbitrary a > 0, this integral is, in general, divergent. As a regularization of the divergent 
integral, the following two alternative definitions for Df exist [Mai96] 


RLD^o,t)fit) = DJtLfit) = 




1 


f{T)dT 


r(n — a) df^ Jq {t — r)“+i 

1 /•* /(")‘i^(T) 


(1.104) 


(1.105) 


r(n — a) Jq {t — t)“+i " 

where n—l<a<n, n = l,2,.... Eq. (11.1041) is the Riemann-Liouville derivative, while 
Eq. p.1051) is the fractional derivative in the Caputo form [Mai96]. Performing integration 
by part in Eq. (I1.104p and then applying Leibniz’s rule for the derivative of an integral and 
repeating this procedure n times, we obtain 


n —1 


DRLfit)=D^f{t) + J2f^''H0+) 


-ik—c 




T{k — a -|-1) 


The Laplace transform can be obtained for Eq. ()1. 10511 . If C[f{t)] = /(s), then 


n—1 


^ mm] = 


-l-k 


k^O 

The following fractional derivatives are helpful for the present analysis 


mm = 


r(i - a) 


, D^[1] = 0. 


(1.106) 


(1.107) 


(1.108) 
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We also note that 


DRLt^ = 


1 ) 

r(/3 + 1 - a) ’ 


(1.109) 


where /3 > — 1 and a > 0. The fractional derivative from an exponential function can be 
simply calculated as well by virtue of the Mittag-Leffler function (see e.g., [Pod98, BH55]): 


E^^iz) = 


oo u 


S) 


Therefore, we have the following expression 


( 1 . 110 ) 


( 1 . 111 ) 
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